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Many ion species commonly used for laser-cooled ion-trapping studies have a low-lying metastable 
state that can become populated due to spontaneous emission from the excited state. 

This requires a repumper laser to maintain the ion in the Doppler cooling cycle. Typically, the 
state, or some of its hyperfine components if the ion has nuclear spin, has a higher multiplicity 
than the upper state of the repumping transition. This can lead to dark states, which have to 
be destabilized by an external magnetic field or by modulating the polarization of the repumper 
laser. We propose using unpolarized, incoherent amplified spontaneous emission (ASE) to drive the 
repumping transition. An ASE source offers several advantages compared to a laser. It prevents 
the buildup of dark states without external polarization modulation even in zero magnetic field, it 
can drive multiple hyperfine transitions simultaneously, and it requires no frequency stabilization. 
These features make it very compact and robust, which is essential for the development of practical, 
transportable optical ion clocks. We construct a theoretical model for the ASE radiation, including 
the possibility of the source being partially polarized. Using ^^Sr"*" as an example, the performance 
of the ASE source compared to a single-mode laser is analyzed by numerically solving the eight- 
level density-matrix equations for the involved energy levels. Finally, a reduced three-level system 
is derived, yielding a simple formula for the excited-state population and scattering rate, which can 
be used to optimize the experimental parameters. The required ASE power spectral density can be 
obtained with current technology. 

PACS numbers: 32.80.Xx, 42.50.Gy, 37.10.Rs 



I. INTRODUCTION 

Optical frequency standards based on single trapped 
ions already have accuracies better than that of the pri- 
mary cesium fountain standards [H-Qi and the Interna- 
tional Committee for Weights and Measures (CIPM) rec- 
ommends quadrupole transitions in three ions (^^Sr"*", 
i99jjg+^ and ^''^Yb"'") as secondary representations of 
the second 01 ■ Amendments to this list have been pro- 
posed by the joint consultative working group on fre- 
quency standards. To characterize the optical clocks 
at the lowest-obtainable uncertainty level, local compar- 
isons between clocks from different institutes are needed. 
This requires at least one of the clocks to be trans- 
portable. Compared to standard laboratory setups, this 
places stringent conditions regarding compactness and 
robustness on all parts of the clock, including the light 
sources. 

Many ions commonly used for ion clocks have an alkali- 
metal-like atomic structure. The analysis in this paper 
focuses on the ^^Sr"*" ion, but the applicability to other 
ions will be discussed in Sec. IVIII The trapped ion is 
Doppler cooled and fluorescence detected using a cooling 
laser driving the ^>S'i/2- transition (see Fig.[T|). The 
^5*1/2- ^1)5/2 quadrupole transition is used as the refer- 
ence or clock transition. The ^Pi/2 excited state has a 



finite probability of decaying to the low-lying metastable 
^1)3/2 state and a repumper laser tuned to the ^1)3/2- 
transition is required to maintain the ion in the 
cooling cycle. As the multiplicity of the "^0^/2 state is 
higher than that of the state, a single-mode re- 

pumper laser will create dark states within the ^^3/2 
sublevels if the external magnetic field is weak Q, 
which is usually required in ion clocks. This causes the 
scattering rate to drop to nearly zero, which inhibits both 
the cooling and the detection of the ion. Dark states can 
be prevented by modulating the polarization of the re- 
pumper [1, Q, e.g., by using an electro-optic modulator 
(EOM). In addition, the laser must be frequency stabi- 
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FIG. 1. Lowest energy levels and relevant transitions for ions 
Ca"*" (n = 4), Sr"*" (n = 5), and Ba"'" (n = 6). Items in gray 
apply to Yb+ only [n — Q). 
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lized. Alternatively, one can use a multimode fiber laser 
and overlap two beams with different polarizations [lo| . 
which results in a fluctuating net polarization. 

We propose using unpolarized, incoherent amplified 
spontaneous emission (ASE) to drive the repumping 
transition. The ASE source offers several advantages 
compared to a laser repumper. It prevents the buildup 
of dark states without polarization modulation using a 
single beam even in a zero-magnetic-field environment. 
Due to the broad bandwidth of the source, it requires no 
frequency stabilization and, in ions with nuclear spin, it 
can excite multiple hyperfine transitions simultaneously. 
These features, and the fact that the ASE source can be 
fiber based, make it very compact and robust, which is 
essential for the development of practical, transportable 
optical clocks. A similar ASE source at the ^-05/2-^^3/2 
transition wavelength could be used to replace the so- 
called clcarout laser, whose purpose is to shorten the 
measurement cycle by emptying the metastable ^1)5/2 
state after a quadrupole transition has occurred. 

In the context of ion traps, incoherent light sources in 
the form of UV light-emitting diodes (LEDs) ha ve p re- 
viously been used for photoionization of calcium |lll4l3l | 
and barium p^ . 

This paper is organized as follows: In Sec. |lll we con- 
sider two well-known models for stochastic fields and con- 
struct a model with adjustable intensity correlations for 
the ASE field. In Scc. lIIIl we solve the eight- level density- 
matrix equations of the ion numerically and compare the 
different field models with each other and with a single- 
mode laser repumper. The effect of a partially polarized 
field is analyzed in Sec. IIVI and in Sec. |V] we derive a 
reduced three-level system that can be solved exactly. In 
Sec. I VII some experimental questions are addressed and in 
Sec. IVIll we discuss which other ion species this concept 
can be applied to. Finally, the results are summarized. 



II. AMPLIFIED SPONTANEOUS EMISSION 
FIELD MODEL 

The electric field of the ASE light at the position of 
the trapped ion can be represented by the sum of two 
components with orthogonal polarizations (ui • U2 = 0): 
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[ui£i(t) + U2f2W]e- 
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where wase is the center frequency of the field and c.c. 
stands for complex conjugate. If the light is unpolarized, 
the components £i(t) and £2it) are completely uncorre- 
lated stochastic variables and their intensities are equal 

We use Ornstein-Uhlenbeck (OU) processes [3 to 
model the stochastic variables, which is convenient as 
there are exact updating formulas for both the OU pro- 
cess and its time integral [l7|. The OU process X{t) 



obeys the Langevin equation 



(2) 



where and c are positive constants, traditionally re- 
ferred to as the relaxation time and diffusion constant, 
respectively, and r(i) is delta-correlated Gaussian white 
noise fulfilling (r(t)r(i')} = S{t - t'). The correlation 
function for X{t) is given by 



{X{t)X{t + T)) 
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Two common types of stochastic fields are the phase- 
diffusion (PD) field and the chaotic field [l^l- The phase- 
diffusion field. 



(4) 



has a constant amplitude and a fiuctuating phase and 
is typically used to describe a single-mode laser of finite 
bandwidth (l9l - [2^ . In this model, the time derivatives 
of the phases, i.e., the frequency fiuctuations, fulfill the 
correlation function ([3]), 



0) = f e-^l^ij,. 



(5) 



with P = and b = ct^. The delta function indicates 
that the two polarization components are uncorrelated. 
The (normalized) first-order correlation function of the 
phase-diffusion field is [T9| 



(r) = 



exp 



1 



■<Pj{t+T)] 
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This corresponds to a lineshapc that is essentially 
Lorcntzian within ±/3 from line center and falls off as 
a Gaussian outside this region. In the limit /3 ^ 5, the 

g-6kl/2^.^.^ cor- 



corrclation function becomes j,, (T) 



&kl/2A. 

responding to a Lorcntzian line with a full width at half 
maximum (FWHM) of b. In this case the OU process 
approximates a Gaussian white noise (Wiener) process. 
As the amplitude of the PD fields is constant, the second- 



order correlation function is 



.(2) 
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We now turn to the chaotic field. 



(r) = 1. 



(7) 



where the real and imaginary parts are independent 
Gaussian stochastic variables. It has both amplitude and 
phase fluctuations and can be used to describe a mul- 
timode laser field with a large number of uncorrelated 
modes Now the electric field fulfiUs Eq. dSj, 

(£,^(t)£j(t + r))=£2e-%l/25,^.J,,, (8) 
where b = 2/Tr and £q = cti/2. This gives the first-order 



correlation function 5ch\j"(''') 



-&|t 



/^(Sio, identical to 
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that of the PD model in the /3 ^ 6 limit. Thus also 
the spectrum of the chaotic field is a Lorentzian with a 
FWHM of b. For a real Gaussian process, the first-order 
correlation function contains all the information about 
the higher-order correlation functions [see, e.g., Eq. (2.1- 
8) in Ref. [l^l] and the second-order correlation function 
becomes g^i^ijiT) = 1 + e^^l'^l(5jj, which shows photon 
bunching. 

The second-order correlation properties of the polar- 
ized ASE from a superlumincscent light-emitting diode 
(SLED) have recently been studied [24j. It was shown 

(2) 

that 3ase(''') depends on the pump current of the SLED: 
at low currents the SLED emits pure chaotic light, 
Sase(O) ~ ^' current increases, the value 

of 

5ase(0) decreases, being 1.3 at the highest current. 
We expect other ASE sources to show a similar depen- 
dence on the pump and have therefore constructed a 
field model with adjustable intensity correlations. Ex- 
pressing the chaotic field as £ch,j{t) = 2~^/^{[£^j(t)]^ 4- 

[f:](t)]2}l/2e«arg[£j(t)-hz£'(t)] ^ | arg £<.h,, («) ^ 

write the ASE field as 

^ASE,,(t) = [^o' cos2 a + |£ch,,(i)P sin2 a] ^'^ e^-'^^--^*). 

(9) 

We have numerically verified that the first-order correla- 
tion function of also this field is g"^Q,-^ij(j) = e~^l'^l/^i5ij. 
The second-order correlation function can be evaluated 
in closed form and is 5AgE ij(''") = 1 + "svv^ aeT^^'^^bij. 
Hence the spectrum is a Lorentzian with a FWHM of 6, 
while the intensity correlations can be adjusted using the 
parameter a. 

Note that while Eq. ^ is identical to the chaotic field 
iS'chj"(^) for a — 7r/2, it will not obey the same statistics 
as the PD field £'pDj (i) when a = 0. It should rather be 
interpreted as a chaotic field with suppressed intensity 
correlations. This is not considered a problem, as the 
ASE source will never fully exhibit a single-mode-laser- 
like behavior. In addition, the repumping efficiency does 
not depend on the actual statistics, as will be shown in 

secnni 

The correlation functions given so far have been for 
the separate polarization components, as indicated by 
the subscript zj. The first-order correlation function 
of the total, unpolarized field ([J) only difi'ers by the 
factor e"'"*^"^"^. However, the second-order correlation 
function for an unpolarized light beam (for a measure- 
ment that does not distinguish the two polarizations) is 
given by 5(2) (^) = [1 + (r)]/2 [H, where gf (T)_ is 
the second-order correlation function of the polarization 
components. This means that the photon bunching is 
reduced by a factor of two for unpolarized light. 

III. EIGHT-LEVEL SIMULATIONS 

We want to demonstrate that the ASE repumper 
prevents the buildup of dark states in the metastable 
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FIG. 2. (Color online) (a) Eight-level system. Fine structure 
state designations are given at the left and magnetic quan- 
tum numbers at the bottom. The sublevels are labeled con- 
secutively and optical dipole transitions are shown as arrows 
together with their corresponding relative transition ampli- 
tudes Cij. The spontaneous decay rates are shown next to 
the wavy arrows, (b) Reduced three-level system obtained 
from the eight-level system in (a) as described in Sec. [V] 

state and between the metastable state and the 
^'5'i/2 ground state. We start by numerically solving 
the density-matrix equations of the eight-level system in 
Fig. Ula) using the formalism described in our previous 
work [7|, with the repumper described by Eq. ([l} and 
either Eq. Q or Eq. ([9|) [the chaotic case is obtained 
from Eq. ^ with a = 7r/2]. All numerical results are for 
88Sr+. 

We are interested in the (near-)zero-magnetic-field sit- 
uation, so we choose the quantization axis along the 
linear polarization of the cooling laser that travels in the 
Uj^ direction. The ASE field is assumed to propagate in 
the same direction and a natural choice of the ASE po- 
larization components is then Ui = u^, and U2 = . For 
the cooling laser we use the (two-level) Rabi frequency 
51c/r = 1 (for ^^Sr^ this corresponds to the intensity 
Ic = 39.8 mWcm~2), which we have shown to give a 
high scattering rate with only minor power broadening 
[3], and the detuning Jc/F = —0.5, which minimizes the 
temperature achievable with Doppler cooling. Here F 
is the spontaneous decay rate of the excited state 

(F/(27r) = 21.6 MHz for 88Sr+). 

As the ASE field is a stochastic variable, this is true 
also for the repumper Rabi frequency. Therefore, we 
use its rms value to describe the strength of the field, 
57r = [(|r2i.(t)p)i]^''2, where { )t refers to a time average. 
For a laser repumper with a linewidth ^ F, it is desirable 
to have a repumper Rabi frequency ilr.iascj_,~ f^c = P 
in order to maximize the scattering rate [7[. The ef- 
fective intensity per bandwidth is thus proportional to 
f^?i==„./r = P- For the ASE source, with a bandwidth 
6 ^ P and a Lorentzian lineshape, to have the same ef- 
fective intensity per bandwidth, we require ^"^/b = F, 
i.e., the Rabi frequency should scale as fir/P = {b/VY/"^. 

Using the two-level Rabi frequency = 
{fi/h)y^2lQ/eoc 0, the peak intensity versus total 
power of a Gaussian beam, Iq = 2Ptot/(7''w^), where w 
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FIG. 3. (Color online) Scattering rate as a function of time 
for (a) a PD repumper, (b) a chaotic repumper, (c) an ASE 

repumper with (/^se iji^) ~ l-'^' ^^'^ i'^) three-level solu- 
tion from Sec. FVl [solid (black) curves]. The curves in panels 
(a-c) are averages over 50 simulations. For comparison, the 
scattering rates for a polarization-modulated laser repumper 
[dotted (blue) curves] and for a laser repumper with station- 
ary polarization [dashed (red) curves] are shown in each panel. 



is the beam waist, and the peak power spectral density 
(PSD) of a Lorentzian hne, SP/Sto = 2Ptot/(7r&), we can 
relate the parameters used in the calculations to the 
power per unit wavelength of the ASE source: 
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For ^**Sr+ this is 
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As it is common to use a bandpass filter at the cooling- 
laser wavelength for fluorescence detection, we define the 
scattering rate as Fgc = AgT{pY7 + pss), where Ag is the 
decay probability to the ground state. 

As the density matrix describes ensemble averages, but 
each numerical simulation corresponds to a particular re- 
alization of the stochastic field, we should average the re- 
sults over a large number of simulations. Figure |3] shows 
the scattering rate for small values of t for a repumper 
with b/T — 100 (bandwidth approximately 2.2 GHz or 
0.0086 nm) and fl./T = (&/r)i/2 = 10 described by 
[Fig. 3(a)] the PD, [Fig. 3(b)] the chaotic, and [Fig. 3(c)] 
the ASE models [solid (black) curves]. The curves are 
averages over 50 simulations, so there are some fluctua- 
tions left that would vanish for a high number of averaged 
simulations. The ASE curves are compared to the result 
for a polarization-modulated repumper laser (modulation 
frequency 0.3r and amplitude tt, see Ref. Q for details) 



[dotted (blue) curves] and for a single-mode repumper 
laser with stationary polarization [dashed (red) curves]. 
In both cases the Rabi frequency was 51r.iascr/r ~ 1, as 
discussed above, and the detuning was (5r.iasor/r = 0.5 in 
order to tune the system away from the coherent popu- 
lation trapping (CPT) resonances that otherwise can oc- 
cur between the ground and metastable states. For the 
single-mode laser with stationary polarization, the scat- 
tering rate decays exponentially at the optical pumping 
rate Fop/F = 0.002, whereas the scattering rate for the 
modulated laser stays constant after a small initial decay 
except for oscillations at the modulation frequency (an 
extended plot up to Tt = 2000, comparing a coherent 
repumper with fixed polarization to one with modulated 
polarization, can be found in Fig. 8 of Ref. [3]). 

As can be seen from Figs. [31[a)-[3t^c), the differences 
between the results for the three field models are negligi- 
ble. Except for the residual fiuctuations, the three curves 
are also very close to the three-level solution that will be 
derived in Sec.|Vl shown in Fig.[3][d) [solid (black) curve]. 

In practice, one is mainly interested in the time- 
averaged quasi-steady-state scattering rate, i.e., when 
the only time dependence is due to the fluctuations. 
For the parameters in Fig. [3l both the PD model and 
the ASE model with different values of 5ASEy (0) (i-^-' 
including the chaotic case) yield a scattering rate of 
(7.8±0.2) X 10*^ s-i for ^^Si+ when averaged over 10 sim- 
ulations, with the variations between single simulations 
being larger than the variation between the averages for 
the different models. The corresponding scattering rate 
for the modulated laser is 8.7 x 10^ s~^. i.e., the scat- 
tering rate for the ASE source is only about 10% lower. 
For a nonmodulated laser repumper, the scattering rate 
(10) drops to 80 s^^, i.e., practically to zero. This value can 
also be estimated as the scattering rate obtained for a 
modulated laser multiplied by the ratio 7rn/rop, where 
7m = 2.3 s~^ is the decay rate of the metastable state 
and Fop ~ 0.002F is the rate of optical pumping into 
this state . 

The effects of a PD field and a chaotic field on a two- 
level atom have been compared by Georges and Lam- 
bropoulos [3] ■ They conclude that the difference is small 
when the saturation parameter S = 20'^^^^!^. /T (T + b) <C 
1. Here, Crms is the root-mean-square transition ampli- 
tude that will be derived in Sec.|Vl For the parameters in 
Fig. [3l S = 1/3, but this system is much more complex 
than a two-level atom: the repumper is driving one leg 
of a A system with Zeeman substructure and the main 
mechanism that reduces the scattering rate is the buildup 
of coherent dark states between the sublevels, not satura- 
tion. Our results indicate that the amount of amplitude 
fluctuations and the exact statistics of the field are not 
relevant for the ASE repumper to work. 

Figure [3] was plotted for a bandwidth of b/T = 100. 
However, the unpolarized repumper principle works over 
a wide range of bandwidths. Figure |4] shows the time- 
averaged quasi-steady-state scattering rate as a function 

of the bandwidth for an ASE repumper with g\g^ ,y (0) = 
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FIG. 4. (Color online) Time- averaged quasi-steady-state scat- 
tering rate as a function of repumper bandwidth b/T for an 

ASE repumper with g^lE.ijW = For 6/r < 1, the re- 
pumper Rabi frequency was fir/r = 1 (black circles), whereas 
for b/r > 1, it was fir/F — y/b/T (black crosses). For the 
smallest bandwidths, the corresponding values for the PD 
model are shown for comparison [gray (red) squares]. The 
circles, crosses, and squares are the mean and the error bars 
the standard deviation of ten simulations. 5r/F = 0.5. 



T5 (black circles and crosses). It shows that, for small 
bandwidths, the polarization fluctuations are too slow to 
efficiently destabilize dark states, so the scattering rate 
is low and the differences between individual simulations 
are large. Also shown in Fig. |4]is the corresponding data 
for the PD model for the smallest bandwidths [gray (red) 
squares]. Here, the scattering rate is slightly lower and 
the simulation-to-simulation variations smaller than for 
the ASE model, presumably due to the fact that the 
phase fluctuates much more smoothly in the PD model 
than in the ASE and the chaotic models. On the other 

('2') 

hand, varying ffASE ijX'^) from 1 to 2 in Eq. (|9]) does not 
affect the results of the ASE model significantly, indicat- 
ing that the amplitude fluctuations are of minor impor- 
tance for the destabilization of dark states. For band- 
widths 6/r > 0.01, there are no signiflcant differences 
between the scattering rates for the different models. 

For 6/F « 0.1 . . . 0.3, dark states are effectively desta- 
bilized while the bandwidth is still narrow, and thus 
the scattering rate is approximately the same as for a 
polarization-modulated laser. As h/T is increased above 
1, with the Rabi frequency scaling as fir/F = (6/F)^/^, 
the RMS fluctuations of the scattering rate (in a sin- 
gle simulation) increase, but at b/T k, 10 the fluctua- 
tions have essentially reached their maximum amplitude. 
More importantly, the time-averaged quasi-steady-state 
scattering rate remains essentially constant for 6/r > 1 
(see Fig. |4]). Thus for 6/r 1 the time-averaged 
quasi-steady-state scattering rate depends only on the 
ratio r2^/6, proportional to the PSD at line center. The 
bandwidths of real ASE sources can be of the order of 



1 ... 10 nm, corresponding to 6/r 10* . . . 10^, but sim- 
ulating this large bandwidths is not practically feasible. 

It should be noted that in a single-ion experiment, the 
fluctuations in the scattering rate due to the stochastic 
ASE fleld will not be detected, since the scattered pho- 
ton detection time constant typically is several orders of 
magnitude greater than the time scale of the fluctuations. 
In fact, the time constant is longer than all time scales of 
the ion system and thus one is essentially measuring the 
ensemble-averaged quasi-steady-state scattering rate. 

IV. PARTIALLY POLARIZED FIELD 

In order for the ASE repumper to be a reliable, practi- 
cal device, it is important that it works even if the output 
light is not completely unpolarized. The ASE field ([T]) 
can be made partially polarized by redefining the electric 
field of the second polarization component U2 as 

f2,Pp(t) =Ci£i(t)+C2£2(t), (12) 

where ci and C2 are complex numbers obeying |ci[^ + 
|c2|^ = 1. The degree of polarization [2^] of the total 
field is then given by P = |ci]. If ci is real, the polarized 
part is linear, if it is purely imaginary, it is circular. 

Figure O shows the time-averaged quasi-stcady-state 
scattering rate as a function of P. The ASE model with 

g\i,^ijifi) = 1-5 and 6/r = 100, and a linear partial 
polarization were used for these simulations, but we have 
found no significant differences for the other field models, 
other values of g^gg other bandwidths, or a circular 
partial polarization. The scattering rate remains within 
10% of the maximum value up to P « 0.8, demonstrating 
that the ASE repumper is not sensitive to the output 
being partially polarized. For fully polarized light, P = 1, 
we obtain the same scattering rate Fgc = 80 s~^ as for 
an unmodulated laser. The inset in Fig. [5] shows the 
scattering rate as a function of the external magnetic field 
for P = 1. We note that if the ASE light is completely 
polarized, the magnetic field required to destabilize dark 
states is of the same magnitude as for a laser repumper 
[3 , regardless of the fact that the ASE is incoherent and 
broadband. 



V. REDUCED THREE-LEVEL SYSTEM 

The numerical eight-level simulations described in 
Sec, mil are very computer-time consuming. We arc there- 
fore motivated to derive a simpler, approximative system 
that can be used, e.g., to find the optimum power spec- 
tral density (PSD) of the ASE repumper for different 
cooling-laser Rabi frequencies and detunings. 

The purpose of the ASE repumper is to provide ef- 
ficient repumping to the cooling cycle while preventing 
coherences from building up between sublevels of the 
metastable state and between sublevels of the ground and 
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FIG. 5. Time-averaged quasi-steady-state scattering rate as 
a function of the degree of polarization for an ASE field with 

5ase iji^) ~ ^^'^ ^/^ ~ 100- ^^^^ circles are the mean and 
the error bars the standard deviation of ten simulations. The 
inset shows the scattering rate as a function of the external 
magnetic field B for P = 1. 



metastable states. The numerical simulations in Sec. IIIII 
show that these coherences fluctuate around zero with 
small amplitudes. Hence they vanish in an ensemble av- 
erage and we set the metastable coherences Pmm' and 
the ground-metastable coherences pgm to zero in order 
to simplify the system. Also justified by the numerical 
results, we set the excited state coherence prs to zero, as 
the ASE repumper is incoherent and the linearly polar- 
ized cooling laser only drives the |1) \7) and |2) —5- |8) 
transitions without coupling the two excited-state levels 
to each other. With the previous simplifications, one no- 
tices that the density matrix elements pis, P27, and pi2 
only depend on each other. As they initially are zero, 
they will remain so and can be discarded. 

The m-e optical coherences decay at the rate h/2 ^ F. 
We can therefore use the "broad-line approximation" [22], 
W\ , which assumes that the optical coherences follow the 
populations adiabatically so that we can use the steady- 
state expressions obtained by setting dpme/dt = 0, 



where Cem is the relative transition amplitude, /* is the 
complex conjugate of the polarization component ampli- 
tude, and q — 0, ±1 is the polarization that connects 
the two |m) and |e) sublevels (see Ref. Q for details). 
Substituting Eq. (|13p into the density-matrix equations, 
we are left with ten equations: the eight populations and 
the two g-e optical coherences pn and p2s- 

Because of symmetry, we must have pn = P22 and 
P77 = P88- As C71 = —Cs2 and the two transi- 
tions |1) — > |7) and |2) |8) are driven by the same 
polarization component, it follows from Eq. (|13p that 
pi7 = —p2s- In addition, from the numerical simu- 
lations we know that the metastable-level populations 
are approximately equal. We can thus define the to- 
tal populations pgg = pn + P22, Pmm = Yfk=3Pkk, 
and Pee = Pn + Pss- If we furthermore formally de- 
fine pge — 2p28 ~ P28 ~ Pi7i wc obtain the equations 
for the reduced three-level system [in the rotating wave 
approximation and rotating frame, see Fig.[2Jb)], 



Pe 



-Tpe 



'-'rms^ 'r / Fmm 



Pmm — -^m^ Pee ^mPmm 
Pgg Pee ^~ ^mPmm 



2 ,^eej - C82r!cPp,4a) 

"-^rms^'r / Hm 



h V 2 

Cs2^cP' 



/Oee()4b) 

(14c) 



Pee), 



(13) 



/F \ C f2 

Pge^ - i-^+iSA Pge+i-^J-^{Pee- Pgg)- (14d) 

Here, pj^p = lm{pge), Ag and A„i are the decay probabili- 
ties to the states \g) and \m), respectively {Ag + Am = 1), 
7,„ is the decay rate of the metastable state, and Sc is 
the cooling laser detuning. The root-mean-square transi- 
tion amplitude Crms is obtained from C^^^ = |/iiC86p + 

mCs5? + \f+M^ = l/ilC75p + |/5C74p + |/;iC73p = 

1/6; where the average polarization probabilities for an 
unpolarized plane wave are l/LiP = = 1/4 and 

l/5f = 1/2. 

Note that Eqs. depend only on the ratio ^"^/b, 
which was already observed for the eight-level results in 
Sec. IIIII These equations can be numerically integrated to 
obtain the time-dependent solution and Fig.[3]shows that 
the solution agrees very well with the numerical eight- 
level results. 

The steady-state solution of Eqs. [M] can be solved ex- 
actly to give the excited state population 



(15) 



where R = C'^^^^VL^/b. Since the metastable state decay further be simplified to 
rate 7^ is much smaller than both F and i?, Eq. (|15p can 



A„ 



f4-h2A 



^1 

■R) 



(16) 
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The last term in the denominator of Eq. (|T6|) describes 
optical pumping and is similar to what one obtains if the 
populations of the \g) and |m) sublevels are equalized by 
a relaxation process such as collisions [2^. Note, how- 
ever, that the exact solution (|15p simplifies to the famil- 
iar expression for a two-level system in the limit Ag ^ 1 
{Am -T" 0) and R ^ 0, which the approximation (fT6|) 
does not. 

For the parameters used in Fig. [3l the three-level 
steady-state scattering rate F^c ~ AgTpee is 8.0 x 10^ s~^, 
in good agreement with the numerical result of (7.8 ± 
0.2) X 10^ s~^ obtained in Sec. IIIII Figure [6] shows the 
scattering rate as a function of both ilf./{br) and the 
experimentally more accessible PSD 6P/6X (assuming a 
beam waist w = 50 yum) for three different cooling laser 
Rabi frequencies. The figure shows that for each value of 
ilc, there is a limiting ASE PSD, above which the scat- 
tering rate does not increase anymore. As anticipated in 
Sec. Illli this value is fl^/b w ft^/T. However, note that 
one comes very close to the maximum scattering rate al- 
ready for PSDs one order of magnitude lower than this 
value. 



n^/br 

r 

,8 lO"* 10"^ 10"" 10"' lo' 




SP/81 (nW/nm) 

FIG. 6. (Color online) Scattering rate as a function of n?/(br) 
(top axis) and PSD (bottom axis, assuming a beam waist 
w = 50 fim) for three different cooling laser Rabi frequencies: 
fic/r = 3 [solid (black) curve], fJc/P = 1 [dashed (red) curve], 
and Slc/r = 0.3 [dotted (blue) curve]. 



VI. EXPERIMENTAL CONSIDERATIONS 

A sufficient PSD at the Sr"*" repumping wavelength of 
1092 nm can be obtained from the ASE of a laser-pumped 
ytterbium-doped fiber [HI. This ASE source consists of 
a pump diode laser and standard fiber optic components. 
The fact that it requires no frequency stabilization or ex- 
ternal polarization modulation not only provides greater 
simplicity and reliability in operation, but also makes it 
a cost-effective alternative to, e.g., a custom-ordered dis- 



tributed feedback (DFB) laser. The emitted light is un- 
polarized and its transversal coherence is similar to that 
of a laser, so it can be focused to the spot size required 
for an ion trap. 

Current SLEDs can deliver a few hundred yuW/nm at 
1092 nm. As these typically emit nearly fully polarized 
light, two SLEDs have to be combined with orthogonal 
polarizations. Taking into account losses and the need for 
filtering (see below), this PSD is at the lower limit of ef- 
ficient ASE repumping. However, more powerful SLEDs 
are expected to become available in the future. 

One important question that needs to be addressed is 
whether there are any drawbacks related to the high to- 
tal intensity of the ASE field. If one uses a bandpass 
filter at the cooling-laser wavelength for fluorescence de- 
tection and ion imaging, these will not be affected by the 
repumper light. The ASE spectrum might have a broad 
pedestal around the peak at 1092 nm. If this is the case, 
light at the clearout wavelength 1033 nm (^£'5/2-^^3/2 
transition) has to be filtered out, otherwise the repumper 
interferes with the detection of clock transitions. 

Due to the high total intensity of the ASE light, the 
light shift of the reference transition is large. If we assume 
an ASE spectral width of 2 nm (FWHM) , corresponding 
to b/r « 2.3 X 10^ and require n,/T = (5/F)i/2 ~ 150, 
the light shift is —14 kHz. The ASE repumper light must 
therefore be carefully blocked during the clock pulse. The 
ASE source can be switched off electronically by turning 
off the pump (laser) current in order to further increase 
the extinction ratio of the mechanical shutter or acousto- 
optic modulator (AOM) used for switching. For the ASE 
source parameters above, the light shifts of the cooling 
and clearout transitions are 4 and 36 kHz, respectively, 
and can be neglected compared to the natural lincwidth 
of these electric dipole transitions. 

The transition strength of the clearout transition (see 
Fig.[T]) is similar to that of the repumping transition, and 
an ASE source could be used also to replace the clearout 
laser that is used to return the ion to the cooling cycle 
after a quadrupole transition has occurred. In this case, 
the main advantage is that no frequency stabilization is 
required. 

VII. APPLICATION TO OTHER IONS 

If we first consider ions without hyperfine structure, 
the scheme described for ^^Sr"*" in Sec.[T]applies also to the 
commonly used alkaline-earth-metal ions ^^Ca"*" [s^l (re- 
pumper wavelength 866 nm ) and i38Ba+ [U (650 nm), 
as well as to other less abundant even isotopes of these 
elements. Also in these ions, the clock transition is the 
^5*1/2- ^I?5/2 quadrupole transition. The Yb^ ion has a 
similar energy level scheme, except that its lowest ex- 
cited state is the extremely metastable ^^7/2 state. For 
Yb^ , the upper state of the repumping transition is usu- 
ally the [3/2] 1/2 state [H (see Fig. [T]) and an ASE 
repumper at 935 nm could be used for ^^^Yb^ both for 
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the ^-^5/2 quadrupole clock transition [s^ and the 

'^Sxj^-'^F'j ji octupole clock transition (3^ . 

Odd isotopes with nuclear spin J, thus featuring hy- 
perfine structure, also have dark states within the ^133/2 
state. In ^''Sr'*' with / = 9/2, there are dark states in 
the F = h and F = % hyperfine levels that require polar- 
ization modulation [sl]. In addition, due to the hyper- 
fine structure, there are three repumping transitions that 
need to be driven. This has been solved by frequency 
modulating the repumping laser over the few hundred 
MHz range of the hyperfine structure [3^. In addition 
to destabilizing the dark states, a broadband ASE re- 
pumper would solve this problem by driving all allowed 
^^3/2" ^^'i/2 hyperfine transitions. 

Of the odd isotopes, ^""^Yb^ is particularly promising 
for optical clocks, as it has a minimal hyperfine struc- 
ture due to its nuclear spin of 1/2 and, in addition, 
has three possible clock transitions. Here, an ASE re- 
pumper could be used with the Sii2-'^D^i2 quadrupole 
and ^S'i/2- ^^7/2 octupole [1,[33| clock transitions. At the 
National Physical Laboratory (NPL), a magnetic field 
around 20 /iT is used to destabilize dark states in the 
^1)3/2 manifold and two separate repumper lasers are 
used to empty its two hyperfine levels F = \ and F = 2 
[13. At Physikalisch-Technische Bundesanstalt (PTE), 
a high laser intensity is used in order to deplete both the 
hyperfine levels and a 1 mT magnetic field, applied dur- 
ing the cooling period, destabilizes dark states in both the 
^£'3/2 manifold and the ^Si/2iF = 1) state [11]. Again, 
a broadband ASE repumper could drive all the allowed 
repumping transitions without creating dark states in ar- 
bitrarily low magnetic fields. An ASE repumper can- 
not be directly utilized with the third clock transition, 
'^Si/2{F = 0) - '^D^/2{F = 2), as the repumper in this 
case must be able to selectively address the ^D^/2{F = 1) 
leve l in order for the electron shelving technique to work 

The oscillator strengths of the repumping transitions 
in Ca"*" [13 and Ba"*" [i^ are comparable to that in Sr+ 
[4l| . so the required PSDs are also similar. The oscillator 



strength of the ^D^/2- ^-D[3/2l i /^ transition in Yb"*" is one 
order of magnitude lower |42l . |43| , and thus the required 
PSD is an order of magnitude higher. Naturally, the 
availability of suitable ASE sources has to be investigated 
for each wavelength. 

VIII. CONCLUSIONS 

We have proposed using unpolarized, incoherent ampli- 
fied spontaneous emission (ASE) to drive the repumping 
transition in a trapped and laser-cooled single ion. We 
have constructed a theoretical model for the ASE radia- 
tion and analyzed its performance compared to a single- 
mode laser repumper by solving the eight-level density- 
matrix equations for a Sr"*" ion. The performance of the 
ASE repumper depends mainly on its power spectral den- 
sity at the repumping wavelength and it works even if the 
field is partly polarized. We have also derived a reduced 
three-level system that can be solved exactly, enabling 
rapid comparison with experiments and optimization of 
experimental parameters. The required ASE power spec- 
tral density can be obtained with current technology and 
the idea can be applied also to the clearout transition 
and to several other commonly used ions. 

An ASE repumper prevents dark states from forming 
without external polarization modulation even in zero 
magnetic field. In addition, it requires no frequency sta- 
bilization due to the broad bandwidth. This makes it 
robust and compact, which is particularly important for 
transportable ion clocks. 
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